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Abstract 

We study fermionic one-matrix, two-matrix and D-dimensional gauge invariant 
matrix models. In all cases we derive loop equations which unambiguously deter- 
mine the large- iV solution. For the one-matrix case the solution is obtained for an 
arbitrary interaction potential and turns out to be equivalent to the one for the 
Hermitean one-matrix model with a logarithmic potential and, therefore, belongs 
to the same universality class. The explicit solutions for the fermionic two- matrix 
and D-dimensional matrix models are obtained at large N (or in the spherical ap- 
proximation) for the quadratic potential. 
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1 Introduction 



Matrix models are usually associated with discretized random surfaces and 2D quantum 
gravity The simplest Hermitean one- matrix model corresponds to pure gravity 
while a chain of Hermitean matrices describes 2D gravity interacting with c < 1 matter 
[|3|. [3], |5| . A natural multi-dimensional extension of this construction has been proposed 
recently || in connection with induced lattice gauge theories. However, the Hermitean 
matrix models possess, as is well-know, the c = 1 barrier above which the stringy phase 
does not exist. For this reason it is interesting to find out what happens for alternative 
matrix models. 

In the present paper we consider fermionic matrix models which involve matrices 
with anticommuting elements. The simplest one — the fermionic one-matrix model — is 
defined by the partition function 

Z = J d^e- NtrVi " m \ (1.1) 

where \I> and \l/ are the N x N matrices whose matrix elements are independent anticom- 
muting Grassmann variables and V stands for a generic even potential 

oo 

y(v^*) = X>(M)*. (1.2) 

k=0 

A D-dimensional extension of the model ( |1 . 1| ) — the adjoint fermion model (AFM) 
— is defined by the partition function j7| 

Zafm = / UdU^Hd^d^e- 8 ^^ (1.3) 

X,fl X 

with x labeling the sites of a D-dimensional lattice. Here Sp[^,^,U] is the lattice fermion 
action 

S F [V,V,U] = ^Nti 

X 

-cf^^p-^^^+^t^ + ^p+^t^)^^^)]) (L4) 
n=i 

with 

J? = r± 7 „ (1.5) 

being the projectors. The case r = corresponds to chiral fermions while r = 1 is 
associated with Wilson fermions. As is well-known, the chiral fermions describe 2 D flavors 
in the naive continuum limit while Wilson fermions are associated with 1 flavor. 

For D < 1, which is associated with the case of a matrix chain, the gauge field U^ix) 
can be absorbed by a (local) gauge transformation of ty x and ^ x so that the model (|1.3|) 
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reduces to fermionic multi-matrix models. The two-matrix model is associated, in par- 
ticular, with D = 1/2. In contrast to their bosonic counterparts, the fermionic matrix 
models are not studied and the proper critical index ^ str ing is not calculated. This is one 
of the motivations for studies of the fermionic matrix models. 

Another motivation for studies of AFM is a recent interest in the problem of inducing 
QCD at large N which is caused by Ref. 0. There are some important features which 
differ AFM from its Hermitean analogue which is known as the Kazakov-Migdal model: 

i) There is no asymptotic freedom of the gauge coupling constant for AFM with chiral 
fermions or Kogut-Susskind fermions, the latter are associated with 4 flavors in 
the continuum. For this reason the kinetic term of the gauge field is not essential 
similarly to the case of quantum electrodynamics || . 

ii) AFM with the quadratic potential has no instability at any values of m. 

iii) AFM with the quadratic potential undergoes ]7j a first order large- N phase transi- 
tion with decreasing m which is associated with the restoration of area law. 

iiii) The matrices \& can not be diagonalized by a gauge transformation so that saddle 
point methods can not be applied straightforwardly at large N. We shall find, 
however, a large- N solution of AFM by means of loop equations. 

In the present paper we derive the set of loop equations which determines the large- N 
solution of AFM with an arbitrary potential at any D. For the fermionic one- matrix model 
we find the solution explicitly and show that it coincides to any order of the 1 / iV-expansion 
with that of the Hermitean one-matrix model with a logarithmic potential. These two 
models belong, therefore, to the same universality class. For the fermionic two-matrix 
model we find explicitly the large- N solution for the quadratic potential ( |1.2| ) and obtain 
an algebraic equation of degree 2J which determines the solution for the polynomial 
potential of degree 2J. We consider the loop equations of the two-matrix model to any 
order of the 1/iV-expansion and discuss their relation to integrable hierarchies. For the 
.D- dimensional AFM we obtain the set of equations which determines the large- N solution 
for an arbitrary potential and write down the explicit solution in the case of the quadratic 
potential. 



2 Fermion one-matrix model 

2.1 Loop equation for fermion one- matrix model 

The simplest example of fermionic matrix models is the one-matrix model defined by the 
partition function (01). The connected correlators ( — tr (#\I') A:i . . . tr (\I>\]/) fcm ) , 

\ iv iv / conn. 

where the average is defined with the same measure as in ( |1 . 1| ) , can be obtained differen- 
tiating log Z w.r.t. g kl , . . . , g km . Note that the correlators with sufficiently large ki, . . . , k m 
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vanish for finite N identically due to nilpotence of Grassmann variables. We introduce 
the generating functions (loop correlators) 

\ 1 ml conn. 

which are understood as the power series in A^ 2 , . . . , A~ 2 (for finite N they are polynomi- 
als). The multi-loop correlators can be obtained from the partition function by application 
of the loop insertion operator 

where 

5 00 d 

mx) = ~E x to' (2 - 3) 



The loop equation can be derived using the invariance of the integral 

A 

a 2 - m> 

under an infinitesimal shift ^ — > \I/ + £ and reads 



cMd® e ~ NtI — — — — = (2.4) 



The contour G\ encircles anticlockwise the singularities of W(u). Eq.( [2.5|) is supplemented 
with the asymptotic condition 

W{\) -> \ as A -f 00, (2.6) 
A 

which is a consequence of the definition Q2.1). Notice that one obtains the single functional 
equation for W(X). This is due to the fact that trV(VWfy) contains a complete set of 
operators. The m-loop correlator can be obtained from W(X) by (m— l)-fold application 
of S/5V(Xi) according to (gg). 

The loop equation (|2.5| ) can be represented as a set of Virasoro constraints imposed on 
the partition function. The coefficients of expansion of Eq. Q2.5|) in 1/A 2 can be rewritten, 
using (|2.2|) and (|2.3|) , as 

L n Z = 0, n > (2.7) 

where the operators 

^ = E^7^ + 4Ejf4 + (2.8) 

^ dg n+k N 2 dg k dg n ^ k og n 

satisfy Virasoro algebra 

[L n , L m ] = (n - m)L n+m . (2.9) 

The Virasoro constraints Q2.7|) look similar to the ones for the complex one-matrix model 
(see Ref. for a review) where the last term on the r.h.s. of Eq. ( |2.8| ) was absent. 
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2.2 Solution at N = oo 



The loop equation Q2.5| ) can be solved order by order in 1/N 2 . The third term on the 
r.h.s. is suppressed by the factor 1/N 2 . Therefore, one can omit it as iV — > oo so that one 
gets 

Doing the integral on the l.h.s. by calculation of the residues at uj = A and at uj = oo, one 
obtains the quadratic equation for W(X) whose solution reads 



1 1 1 /l ^ 2 



W(X) = -V*(A) + (jAV'(A) + lj " A 2 Q(A), (2.11) 

where Q(X) is an even polynomial of degree 2 J — 2 if V(A) is the one of degree 2 J. The 
coefficients of Q(X) are unambiguously determined by imposing the analytic structure of 
W(X) and by the asymptotic condition ( |2.6| ). 

The explicit solution for the simplest Gaussian potential is 

W(X) = l - Ql X + i - yi g *\* + l. (2.12) 

This solution as a function of A 2 has one cut along imaginary axis which represents the 
'support of the eigenvalue density' of the matrix tyty. Note that this cut does not lie on 
the real axis. As it is shown below, this is the general property of the one-cut solution. 
Another peculiar property of the fermionic Gaussian model is that ( S (fyfy) 2k ) = to 
the leading order of the 1/iV-expansion. The fact that ( || \£ r \E'\E , \E') vanishes can be easily 
verified using, for example, the Wick theorem which results in two terms with opposite 
signs. 

The general solution ( |2.11| ) has n cuts on the complex A 2 plane with n < J for 
V(A) being the polynomial of degree 2J. The one-cut solution can be represented in the 
following form 



W(\) 



du \V'{u) + l 



(* -*)(*-„) 



ICi 47rz A — uj \ {uj 2 — x)(uj 2 — y) ' 
The ends of the cut are determined by the asymptotic condition Q2.6|) which yields 

duj \V\uj) + 1 _ q 

Ol 47TZ ^(^2 _ x )(^2 _ y ) 

duj \uj 2 V'{uj) + 2uj = i 



The solution ( |2.13| ), ( ^.14| ) reminds the simplest one-cut solution (as a function of A 



2^ 



of the Hermitean one-matrix model. This is not a coincidence. As is shown in Appendix A, 
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the fermionic one-matrix model ( |1.3|) is equivalent to all orders of the 1/iV-expansion to 
the Hermitean one-matrix model with the logarithmic potential 

V H ($) = V(V$) + 21og$ (2.15) 

where $ = is the N x N Hermitean matrix. The large- N solution which is associated 
with the fermionic model is, however, somewhat different as compared with the standard 
solutions of the Hermitean one-matrix models, because it is known that the integer value 
of coupling constant in front of the logarithm is critical JT(| for the large- iV limit. In 
particular, the points x and y cannot lie on the real axis. Otherwise the solution were 
have a pole at zero, so that the contour C\ would encircle the origin. One can readily 
verify that Eqs. ( |2.14| ) would have no solution in this case. 

The multi-loop correlators can be calculated using the equation (|2.2|) . For the two- loop 
correlator we have 



W(v, A) 



N 2 A(\ 2 -u 2 ) 2 



(A 2 - x){u 2 - y) + (A 2 - y)(y 2 - x) 



(2.16) 



>(\ 2 -x)(\ 2 -y)(v 2 -x)(v 2 -y) 

which looks like the two-loop correlator [|11| of the Hermitean one-matrix model. 

Eqs. ( |2.14j ) simplify for reduced potentials with g 2 k = which change sign under the 
chiral transformation tyty — > — ^ty. Then for y = —x the first equation reads 

(2k + 1) 



w (X) 



W[X) 



£( 

k 



92k+V 



2 2k (k\y 



-X 



2k+l 



(2.17) 



while the second one turns into identity. Eq. ( 2.17 ) coincides with the proper one for the 
Hermitean one-matrix model with the polynomial potential 



k 



The critical behavior emerges at the point where 



w (x r 



0, 



W(X r 



2 . 



(2.18) 



(2.19) 



This is nothing but the standard multi-critical point of the Hermitean one-matrix model. 
For this reason the critical index jstring coincides in genus zero with the one of 2D gravity: 

1 

~2 ' 



"~) string 



(2.20) 



3 Fermion two-matrix model 



3.1 Derivation of loop equations 

Let us define the partition function of the fermion two-matrix model by 

r _ _ 
Z 2 = / d^ 1 d^ 1 d^ 2 d^ 2 e 



Ntr (-V(\/*i*i)-V(\/*2*2)+c[-*i*2+*2*i 



(3.1) 
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The asymmetric case is associated with non-equal potentials 



and 



(3.2) 



(3.3) 



Let us define the odd-odd and even-even two-point correlators, respectively, by 



G(M) = /£(» 



1 



N 



and 



y (A 2 - * 2 * 2 ) 

A 



(*( V - 

\ N V 2 - Mi) 



(A 2 -^ 2 ) 

Analogously, the correlators of arbitrary powers of ^i^i are determined by 

v 



N Vz/2 _ ^r 1 ^ 1 

which enters the asymptotic expansion of W{v, A) in 1/A: 



(3.4) 



(3.5) 



(3.6) 



n=0 



x ' W 2n (u) 

A 2n+1 



Similarly, we define 



n=0 



G 2n+ i(l/) 



(3.7) 



(3i 



The first loop equation results from the invariance of the measure in 

= 



/ tr / A 1 1 T 

( — (t A = = ^< 

\ N V (z/2 - # # ) - tf 2 tf 2 ) ' 



under the shift 



(*0„+ 



A7 



(3.9) 



(3.10) 



where t A (A = 1, . . . , A^ 2 ) are the generators of the U(N) which are normalized by 



^[t A ]ij[t A ]ki = N5u5 kj 

A=l 



(3.11) 



Quite similarly to Ref. |12j which deals with the Hermitean two-matrix model, this equa- 
tion takes at N = oo, when factorization holds, the form 



duo V(u) 
Ci 47ri v — uj 



G{u, A) = W {v)G(v, A) + c \W(v, A) - W (v) 



(3.12) 
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where the contour C% encircles anticlockwise the singularities of G(u, A) as a function of 



UJ. 



Since Eq. ( |3.12| ) expresses G via W, one needs one more equation which relates these 



quantities. It can be derived making the shift 



(*o«^(*i)«+^(n„ (3-i3) 



A7 



m 

(w^ v-fr.w-W )" - (3 - i4) 

The resulting equation reads 

r ^!M^( W) ^ . ^(^(i/, A) - 2^(i/, A) - dG(i/, A) . (3.15) 

JCl 4:711 V — UJ V 

3.2 Exact solution at N = oo 

Quite similarly to the scalar case |13|], Eqs. ( |3.12| ), ( |3.15| ) can be rewritten as an equation 



for Wq(u). To this aim let us express W 2n (u) and G 2n+ i(i>), which are defined by Eqs. (|3.7|) 
and Op), for n > 1 via W (u): 

-cG 2n+1 {v) = I p,^lw 2n (u) + (-- W (yj) W 2n {y) , 
JCi 4m v — uj \v J 

cW 2n+2 (u) = I ^y^LG 2n+l {uo) - W (u)G 2n+1 (u) . (3.16) 
JCi 4m v — uj 

These recurrence relations are obtained by expanding Eqs. ( p.!2|) , ( p,15|) in 1/A. 



The equation for Wo can now be obtained taking the \ jv term of the expansion of 
Eq. ( |3.12|) (with 1 and 2 interchanged) in \ jv and reads 

J2 ki kG 2k - 1 (X) = c[XW (X) - 1] . (3.17) 

k 

The equations which come from the next terms of the expansion of Eqs. ( |3.12j ) and ( |3.15|) 
in 1/v should be automatically satisfied as a consequence of Eq. ( |3.17|) . For the potential 
V being the polynomial of the highest power 2 J, Eq. ( |3 . 1 7| ) is a 2J-ic algebraic equation 



for Wo (A) similarly to the scalar case \T2, [13 



It is instructive to consider in detail the case of the quadratic symmetric potentials 
V{uj) = V(uj) = t\uj 2 (quadratic potential) (3.18) 

when 

- cG x {\) = h(\W (\) - 1) + (- - Wb(A)) W (X) (3.19) 
and Eq. ( p. 17j ) is quadratic: 

tiGi(A) = c[AW (A) - 1] • (3.20) 
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The solution to ( |3.19|) , ( 3.20|) has the same form as the one ([2.12 ) for the one-matrix 
model with 



9i = h + 



(3.21) 



3.3 Relation to integrable hierarchies 



While Eqs. ( |3.12| ), ( |3.15|) have been derived at N = oo, they can be extended for the 
two- matrix model to any order of the 1/iV-expansion. The difference between N = oo 
and finite- N equations resides in irreducible correlators (of order 1/N 2 ) which can be 
conveniently expressed in the asymmetric case V{v) ^ V(z/) via the loop insertion operator 
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5V(u) 



E 

k=Q UL k 



d 1 



(3.22) 



The set of the loop equations which extends Eqs. (|3.12|) , (|3.15|) to finite N reads 
dU V{UJ) G(u, A) = Wo(y)G{v, A) + JL^—G(u, A) + c[\W(u, A) - W (u)} (3.23) 



iCi Airi v — uj 
and 

duj V{u) 

Ci 47T2 V — UJ 



L 



W(u, A) = W Q {v)W(v, A) + -^tt^W(is, A) - 2-W(v, A) - c\G(v, A) 



N 2 5V{y) 



(3.24) 



so that for the recurrence relations one gets 



-cG 2n+1 {u) = [ ^^W 2n (uj) + (-- W Q (u) ) W 2n {v) - , 
Jci 47T2 v — uj \v J JS^oVyv) 



1 5 



Wo n M 



cW 2n+2 (lS 



duj V{uj) 
Ci 4:7ri v — uj 



G 2n+1 (uj) - Wo(u)G 2n+1 {u) 



I 5_ 



■G 2n+1 {v) .(3.25) 



Introducing the analogues of the loop insertion operator (|3.22| ) for G 2n+ i(v) and 



WoJv) with n > 0: 



£-2n+l( ly 



G 2n +\{y) — —C 2n+ i(u)Z 2 

^2 

OO 1 

,(2n+2) 1 



E w4 2 



k=—n 



V 



2(fc+n+l) 



(3.26) 



and 



W 2n {v) = ^-£ 2n (v)Z 2 , 

^2 

oo -| 

£ 2 „H ^ - £ <" +1, ^ht 

k=—n 



(3.27) 
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so that 

it can be shown after a little algebra that Eq. ( |3.16| ) is now identically satisfied for any 
Z 2 provided that the operators obey the recurrence relations 

oo i k+n Q 

_ cW (- +2 ) = E mtmW ^ + _L E ±-Wt: 1] + k > -n 

m=l JV m=0 ° lm 

cW ( 2 n + i) = ^ mtmW ^) k + _L ^ * W g) • * > -n (3-29) 

m=l m=0 m 

Analogous operators for the Hermitean two-matrix model were advocated by Marshakov 
et al. [HI. 



The set of equations similar to Eqs. ( 3.23 ) to ( |3.2S ) can be obtained by varying the 



partition function (|3.1| ) w.r.t. \&2 rather than ^j. One gets 

1 00 1 

G 2 n+l(A) = - — J2 W k U X 2(k+n+l) Z2 ' 

W^(A) = -- E ^ n+1) ^)TT^ > (3-30) 

with } given by the same formulas as ^ with t m replaced by t m . 

The equation for PVo(A) can be obtained taking the \jv term of the expansion of 



Eq. ( |3.23| ) in \ jv and reads 

E ^ m G 2ro _!(A) = c(A^ (A) - 1) . (3.31) 



m>l 

Similarly, one gets 



E mt m G 2m ^{\) = c(XW (X) - 1) (3.32) 



which determines Wq(\) versus {tk} and {tk}. Eq. (|3.32|) coincides with Eq. (|3.17|) above. 



If all t m 's vanish except for some m = n (t n = l/n), Eq. (|3.32|) reduces to the 
constraints 

ytf n) Z 2 = cWl%Z 2 k>l — n (3.33) 

imposed on the partition function Z 2 . 



4 D- dimensional AFM at large N 
4.1 Arbitrary potential 

Let us consider AFM on a D-dimensional lattice with is defined by the partition func- 
tion ( |1.3j ). Let us define again the odd-odd and even-even one-link correlators, respectively, 
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by 

C#'(i/, A) = ( Sfe o \ r UJx)— J— H+aUi(x)) ) , (4.1) 

where i, j = 1, . . . , s are spinor indices, and 

w ^ - ( I (i^W ^-w^ w)) ' (42) 

Analogously, the correlators of arbitrary powers of ^ x ^x at the same site x are determined 
by 

The asymptotic expansion of W{v, A) in 1/A is defined again by Eq. (|3.7|). 
The first loop equation results from the invariance of the measure in 



x)) ) = (4.4) 



under the shift 

Quite similarly to Ref. [IB] which deals with the Kazakov-Migdal model, this equation 
takes at iV = oo, when factorization holds, the form of loop equation ( p.l2| ) for the 
two-matrix AFM: 

j du^Viul = A) _ l_ / T yi r w( A) _ (46) 

JCi 47U 1/ — LO ^ p S v /L J 

with the potential V given by 

V'(w) = V"(w) - 2m(2D - l)F(u) (4.7) 

where 

a = r 2 -l. (4.8) 

The function 



oo 



F(u) = J2 F nU 2n ~ l (4-9) 

n=l 

in Eq. (|4.7| ) is the one which appears in the one-link correlator 

/ l l(t a U X U j )) =L^A = \ - ■ (4.10) 

\N\ 7 /o.i. r dxdxdU e _JVtr ( v (^)-* p M t/ x^+xP+m^c/]) 

where the averaging is only w.r.t. U and x while if) plays the role of an external field. 
Similarly to the case of the Hermitean model JM , [15| , the following formula holds at 

= oo: 

| (t^xW)) = ( P *T £ (iVW)- 1 ) , (4.H) 

/ o.l. n=1 IN 
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This formula is discussed in Appendix B. 

One more equation can be derived making the shift 

(*0 B -(*0,, 

in 

It takes the form of Eq. ( |3.15| ): 



A7 



(4.12) 
(4.13) 



duo V'H 



A) = WoW% A) - (a + A) - cA (p±) J * G$>, A) . (4.14) 



d 47T2 ^ — UJ 

To solve the system (|4.6|) and ( [4. 14 ), it is convenient to introduce 



G(u,\)=lP±fG«(v,\) 



and 



n=0 A 



Multiplying Eq. (|4.6| ) by (-P^)- 7 * one rewrites Eqs. (4.6) and ( 4.14j) as 

<L ' V UJ) A) = W (u)G(v, A) - ca\\W(y, A) - W„(f) 



and 



ICi Am' v — uj 
duj V(u 



1 



Ww, A) = W (u)W(u, A) - (s + 1)-W(z/, A) - cAG(z/, A) . 

Ci 47T2 z/ — UJ V 



At s = 1 and o = — 1 these equations coincide with (|3.12|) and ( [3.15| ). 
The analogue of Eq. ( |3.16| ) now reads 

-cG 2n+1 (u) = I ^^lw 2n (uj) + - W (vj) W 2n (u) , 



-caW 2n+2 {v) 
while that of Eq. (|3.17|) is 



Ci 4m v — uj v ^ 

rfcu V'(w) 



Ci 47U i/ — uj 



G 2n+1 (uj) - W (v)G 2n+1 (p) 



kt k G 2k -i(\) = -ca[XW (X) - 1] . 



(4.15) 
(4.16) 



(4.17) 



(4.18) 



(4.19) 



(4.20) 



Notice that Eqs. ( f£T7| ) and ( g^g ) (or ( |£ISD and ([00]) ) express A) and G(i/, A) 
via V'(c<j) (or Wo(^)) which is considered as given. To find it for given V'{uo) and, therefore, 
to determine F(u), one can use the equation 



JCi 2m v — uj 
12 



(4.21) 



which results from the definitions (fO]), ( |4.15|) , (|4.16|) and Eq. ( |4.11|) . Using ( |4.19|) 
Eq. ( [4.21| ) can be rewritten as 



du V'(uj) 



Ci 47ri v 



Wn(uj) + 



UJ 



s + l 



<JC 



du F{uj) 

Ci 2lXl V — UJ 



W {uj) 



(4.22) 



which looks similar to the large-iV loop equation (|2.10|) for the fermionic one-matrix 
model. Thus, one concludes that 



2 Cont„ W (u) = V'{v) - 2acF(v) + 



s + l 



(4.23) 



at the cut (or cuts) of W . Eqs. Q477D , ( ^20|) and ( ggg) unambiguously fix V(uj), F(uj) 
and Wo(uj) in the full analogy Jl3| with the Kazakov-Migdal model. For Kogut-Susskind 
fermions when s = 1 and a = — 1, these equations coincide with those of the fermionic 
two-matrix model discussed in the previous section. 



4.2 Explicit solution for quadratic potential 



The analysis for the case of the quadratic potential V(uj) = muj 2 is similar to that of 
Sect. |]| The solution to Eqs. (fl20| ) and (^23|) reads 



Wo(X) = - 



fJ,X + ^ - \j^ 2 \ A + 2(s - + (s + l) 2 

A A 



with 



Analogously, one gets 



(D - l)m + D^m? - Aac 2 (2D - 1) 
{2D - 1) 

2c 



or substituting Eq. (|4.25| ) 



F(uj) 



V/i 2 + 4ac 2 + n 



2c 



UJ 



m 2 - Aac 2 {2D - 1) + m 



■ uj . 



(4.24) 

(4.25) 
(4.26) 

(4.27) 



Eq. ( f4.27p coincides with the result of Ref. [0] where the same formula was obtained 
utilizing the analogy between adjoint and fundamental fermions in the phase with local 
confinement for the quadratic potential. Similarly, picking up the 1/A 3 term in Eq. ( f4.24| ), 

one gets 

(4.28 



N / 



where fi is given by Eq. (|4.25| ), which agrees with the result of Ref. JT6J for lattice QCD 
with fundamental fermions at vanishing plaquette term. Since ( |4.28|) is nonvanishing at 
m = 0, chiral symmetry is spontaneously broken. 
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The fact that ^ || tyty) coincides for adjoint and fundamental fermions is based on 
the sum-over-path representation which is discussed in Appendix B. This coincidence is 
no longer valid for higher correlators, say ^ ^ {^'9) k \ with k > 2 does not coincide with 
the proper correlator for the fundamental representation. The value of this correlator for 
AFM can be obtained picking up the 1/A 2fc+1 of the expansion of (4.24). 



5 Conclusions 

The main result of this paper is that the fermionic matrix models can be solved by the 
method of loop equations quite similarly to the Hermitean ones despite the fact that 
Grassmann matrices can not be diagonalized. We have obtained the explicit solution 
for case of the simplest quadratic potential in the phase with local confinement. It is 
doubtful, however, that a complete set of equations can be written down for only G(u, A) 
and W(v, A), which are defined by Eqs. (|4.1|) and (42), in the phase with area law. 



The fact that the critical index ^string of the fermion one-matrix model coincides in 
genus zero with the one of the Hermitean one-matrix model does not necessarily mean 



that the two models are identical in the continuum. An example is the model of Ref. [fL7 
which is associated with D = open strings. In contrast to the standard case, one might 
expect alternating signs of the genus expansion of the partition function of the fermionic 
one-matrix model which makes it convergent, since the Grassmann integral in ( |1 . 1| ) is 
always convergent. It is hardly to imagine, however, that the double scaling limit of the 
fermionic one-matrix model differs from that of Ref. || by something except signs. This 
continuum limit of the fermionic one-matrix model deserves future investigations. 
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Appendix A Fermion versus Hermitean one-matrix 

models 



We demonstrate in this appendix the equivalence of the fermion one-matrix model and 
the Hermitean one with a logarithm added to the potential. We start from the derivation 
of the loop equations for the Hermitean one-matrix model with the logarithmic potential 



7H 



Me -NtrU(<t>) 



where 



[/($) = y H ($) + a log $ , K H ($) = ]T g k $ k 



(A.l) 



(A.2) 



k=0 



and $ is N x iV Hermitean. The models of this type were used to calculate the virtual 
Euler characteristics of the moduli space of Riemann surfaces [TJ| and are known to 
describe c = 1 matter compactificated on a circle of the self-dual radius interacting with 
2D gravity [H5|. The large- N limit of these models was studied in Ref. \IU . 



To avoid the divergence at $ = 0, let us consider the loop equation which results from 
an arbitrary shift of $ in 



where the average is with the same measure as in (|A.1|) . Introducing the loop correlator 



tr 



one writes the loop equation as follows 



dui wl/' H (wJ 

Ci 2-7T2 A — UJ 



N A - $ 



A(V H (A)) 2 - aW n (X) + 1 



N 2 SV H (\) 



(A.4) 



H/ H (A), (A.5) 



where the contour C\ encircles all singularities of W^(X) and 



5V ll (X) 



k=0 



-i d 



dg k 



is the loop insertion operator. 

Eq. ( |A.5| ) coincides for a = 2 with the loop equation 
model providing 

AH/ H (A 2 ) = W(A) 
where W^(X 2 ) and W(A) are defined by Eq. (|A.4j) and Eq. 



(A.6) 



of the fermion one-matrix 



(A.7) 



1|), respectively, while 



(A.8) 



This proves the equivalence of the fermion one-matrix model and the Hermitean one with 
the logarithmic potential to any order in 1/JV. 
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This equivalence can be formally shown at the level of the partition functions (|1.1|) 
and ( A.l|) . To this aim, let us rewrite ( |1.1| ) asi 



^I^Me-^^'^-W) <x Jd$e- NtlV ^(det$) ™ (A.9) 

which coincides with ( |A.1| ) providing a = 2 and Eq. (|A.8|) holds. The matrix $ is to be 
identified with The integral ( |A.1| ) is, however, ill-defined for finite A" because of the 
divergence at zero. 

The equivalence discussed in this appendix makes it possible to understand why the 
explicit N = oo formulas of Subsect. [T2] look quite similar to the proper formulas for the 
Hermitean one-matrix model (see, e.g., Ref. In particular, the chiral-noninvariant 
case when V involves only odd powers of (^^z) is associated with the Hermitean one- 
matrix model with an odd potential V^(X) = — V H (— A). 



Appendix B Contraction of the one-link correlator 

Let us consider the following correlator 

Gj(0) = ^ (OU^x)^(x + ii)Ul(x))y (B.l) 

where O does not depend on U^(x), Ut(x), ^f(x+fi) and ^/(x + li), but is arbitrary in other 
respects. The averaging is determined by the partition function ( |1.3|) . We shall show that 
the calculation of this correlator within the large mass expansion can be reduced in the 
large- A" limit to the calculation of the correlator which is independent of U^(x), U^x), 
ty(x + ll) and + ll): 



tr / , x\ .-/til FUNf(x)V(x)) » , 

- (O^l^jW^ . [O l ; Y //V (x) ), (B.2) 



X - " ^ ' " > • \ N \ .^( x )^( x ) 

where F(u) is some universal (i.e. C-independent) odd analytic function 

oo 

F(u) = £ F n u 2n ~ l . (B.3) 

n=l 

If we expand the exponent in (|1.3| ) in the power series in all couplings except the 



mass of the fermion field then the calculation of ( |B.1|) reduces to the Wick contractions 
among ^f(y) and ^(y) and group integration over all U^y). The crucial point is that 
the group integration goes independently on each link due to the absence of the kinetic 



term of the gauge field. So it is sufficient to prove ( |B.2| ) for the one-link correlator ( [4.10 ) 



2 Notice that if \P and \& were just complex rather than Grassmann matrices, the determinants would 
cancel while <£> would be positive definite. 
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t a : d x l 

ij> p • : d x' c(p-) p ' 

V s q : » r c(p+)» 



Figure 1 : The graphic representation of the large mass expansion of the one-link corre- 
lator ( gig ). 



The remaining part of the large mass expansion decouples because of the properties of 
U (N) group integrals and the large- N factorization. Although the coefficients F n will be 
different for fl4.10|) and for the full lattice correlator ([B.l] ). 

It is easy to show that ( p.2| ) holds for fl4.10| ). The graphic representation of the large 
mass expansion of ( f4.10| ) for the Gaussian potential is depicted in Fig. [I]. The integration 
over Xi X can De done using the Wick rules (see Fig. after which the group integration 
is trivial, because all U matrices cancel with W. The answer has the form (|B.2| ) with 
F(u>) = F\uj due to the properties of projectors. The value of F\ for the D-dimensional 
model is given by (]427|) . For the non-Gaussian potential the diagrams of the type depicted 
on Fig. |3] appear, but their contributions are also of the form ( |B.2|) , (|B.3|) , i.e. the terms 



like {t°"^)^)^)^)) do not arise. This is because the potential depends only on ipifj. 

In the case of Kazakov-Migdal model the formula ( |B.2|) holds even if the averaging in 
( CT ) is only w.r.t. U^x) [[ 
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This result is nonperturbative (it follows simply from 
the gauge invariance and the large- iV factorization). Also it is possible to calculate the 
pair correlator of the gauge field provided that all the correlators of the type ( ^ $ n ) are 
known [13 . This is not the case for the model under consideration. The reason is that 
one needs to compute a more complete set of correlators than ( || (^ r ^ r ) n ), for example 



the correlator of the type ( | (**^)). 

Of course, our proof of Eq. (|4.11|) is valid only in the local confinement phase which 
is characterized by vanishing expectation values of nontrivial adjoint Wilson loops 

tr 



W A (C) 



N 



U(C) 



5, 



0)^min (C) ? 



(B.4) 



where A min {C) is area of the minimal surface spanned by C. Eq. ( |B.4| ) is valid to all 
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Figure 2: An example of the use of the Wick rules for integration over x 5 X i n t ne 
Gaussian model. 




(P+)^trT^ 



Figure 3: An example of the contraction of color indices in the theory with quartic 
interaction: g is the coupling constant. 



18 



3 



x 



X + fl 



+)J,,X 



Figure 4: The typical paths T x+ ^ jX which contribute to the sum on the r.h.s. of Eq. ( |B.5|) 
in the local confinement phase. These T x+fl)X coincide with L x ^ x+I1 passed 
backward modulo backtrackings which form a ID tree. The sum over T x+fJj ^ x 
is reduced to summing over the backtrackings. 

orders of the large mass expansion and was implicitly used in the proof. 

To demonstrate it, let us consider the following example. The well-known sum-over- 
path representation of the correlator 

(p( X )U»(xM X + riUl( X )) = - £ WA ^Ci\T x) SpP-flW), (B.5) 

^ X-\- f_L,X 

where L x ^ x+fl is the link connecting x with x + and the sum goes over all paths T x+ti x 
complementing L x ^ x+ ^ for a closed contour, P~(T X+ ^ )X ) being the product of ordered 
along this path. In the phase with local confinement one substitutes ( |B.4| ) in this for- 
mula so that the typical paths r x+AtiX coincide with L x+fl)X modulo backtrackings as 
is depicted in Fig. |j. In the phase with normal area law, nontrivial loops contribute to 
the sum in ( p.5|) which were disregarded in our consideration of the group integrals. It 
is doubtful that the formula ( |B.2|) , which enables one to write down a complete set of 
Schwinger-Dyson equations for local objects, really works in the area law phase where 
the dynamics of some extended objects is expected to be nontrivial. 
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